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Abstract Flowing-afterglow measurements in decaying H 3 + or D 3 + plasmas suggest that de- 
ionization does not occur by simple binary recombination of a single ion species. We find that 
vibrational excitation of the ions fails to provide an explanation for the effect contrary to an 
earlier suggestion. Instead, we suggest that collisional stabilization of H 3 ** Rydberg molecules 
by ambient electrons introduces an additional dependence on electron density. The proposed 
mechanism would permit plasma de-ionization to occur without the need for dissociative 
recombination by the mechanism of potential-surface crossings. 

INTRODUCTION 

It is usually taken for granted that typical plasma afterglow experiments 
yield the binary dissociative recombination coefficient for molecular ions and that, 
in principle, the measured values would be the same as those inferred from beam 
experiments. The assumption seems plausible since the capture of an electron into 
a repulsive state of the neutral is thought to induce rapid dissociation 1 , leaving 
little time for third-body interactions. The argument, however, presupposes the 
existence of suitable dissociation channels. If such do not exist, then perhaps long- 
lived Rydberg molecules can be formed and third-body interactions may be more 
important. The experiments described here suggest that this situation may exist in 
plasmas containing H3 + and D3 + ions. 

Theoretical arguments 2,3 indicate that the recombination 

H3+ + e- -> H+H+H 

H 2 +H ( 1 ) 

for vibrationally cold H3 + (v=0) ions by direct dissociative recombination should 
be slow due to a lack of suitable crossings between the ionic and neutral potential 
surfaces. However, vibrationally excited H3 + ions in v >3 are predicted to exhibit 
faster recombination. This expectation disagreed with the early microwave- 
afterglow experiments 4 > 5 which gave fairly "normal" recombination coefficients 
(~ 2 xl 0 ‘ 7 cm 3 /s) for H3 + ions believed to be in v= 0 . This "H3 + dilemma" has 
existed for about 20 years. 

The structure and energy levels of H3 + and H3* are quite well known. A 
brief introduction and extensive references can be found in the "thumbnail sketch 



of H3 + " by Oka and Jagod 6. The importance ot H3 + in astrophysical environ- 
ments is discussed in reviews by Dalgamo 7 and by Tennyson et al. 8 . A fairly 
complete account of the history of dissociative recombination work has been given 
by Bates et al. 9 . 

One possible resolution of the H3 + dilemma was offered by Adams et al. 10 . 
They interpreted their observation of an anomalous decay of a flowing afterglow 
plasma by assuming the presence of two species of H3 + ions, one of which recom- 
bined fairly fast while the other recombined slowly (a<2xl0* 8 cm 3 /s). By identify- 
ing the fast and slowly recombining ions as, respectively, v>3 or v=0 ions, theory 
and experiment could be reconciled. This solution, while being in conflict with the 
microwave afterglow results, seemed quite plausible until Amano 11 ’ 12 showed 
that spectroscopically identified H3 + ions in v=0 recombined readily with electrons 
(a~ 1.8 x 10* 7 cm 3 /s). A repetition of the flowing-afterglow measurements by 
Canosa et al. 13 - 14 gave no evidence for the existence of two distinct types of H3 + 
ions. 

Recombination cross sections for H3 + have been measured repeatedly in 
the merged-beam apparatus at the University of Western Ontario 13, ^ 17, 18 The 
results showed that the cross section varied with conditions in the ion source and 
this was ascribed to vibrational excitation. In later work 1 **, the H3 + vibrational state 
was deduced from the threshold energy for electron-ion dissociative excitation. 
The result indicated that v=0 ions recombined only slowly with electrons. 

When ion storage rings became available, they were very quickly applied 
to H3 + recombination. In these instruments, vibrationally excited ions have time to 
relax by radiative transitions (including the vi=l metastable level of H3 + with a 
lifetime 19 of 1.2 s). Larsson et al. 20 and Sundstrom et al. 21 , using the CRYRING 
storage ring, obtained recombination cross sections for H3 + (v=0) over the wide 
energy range from 0.0025 to 30 eV. The analysis of the experimental cross 
sections yielded a recombination coefficient of (1.15 ± 0.13) xlO' 7 cm 3 /s at an 
electron temperature of 300 K. 

Since the storage-ring measurements also yielded fairly large recombina- 
tion coefficients for v=0 ions, the problem of finding a recombination mechanism 
for v=0 ions became more acute. New mechanisms were proposed 22, 23 that do 
not require crossings between ionic and neutral potential curves. It is possible that 
such mechanisms may play a part in the recombination of H3 + ions, but no 
detailed calculations have been made. 

Not all experimentalists were willing to accept a sizable recombination 
coefficient for H3+ in v=0. Smith and Spanel 24 repeated and refined the afterglow 
experiment of Adams et al.s' 10 and concluded that only H3 + ions in v>3 
reeombined efficiently with electrons. Again, those authors found that the decay 
of H3 + plasmas was anomalous in the sense that the recombination coefficient 
appeared to become smaller in the late afterglow. It was, in part, this work that 
motivated us to carry out an additional experiment. As will be seen, we largely 
reproduced Smith and Spanel's 24 results that there is something peculiar about the 
decay of H3 + plasmas but our interpretation of the experiment is quite different! 
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To avoid ambiguities, we will use the term "recombination coefficient" (a) 
only for purely binary recombination. We will use the word "de-ionization coeffi- 
cient" (P) to denote the quantity that is derived from plasma decay experiments. 
The two quantities may be equal, but this is not assumed from the onset. 

A more complete description of the experiment, additional details of the 
analysis, and the interpretation will be given in a journal article 25 . 

EXPERIMENTAL METHODS 

Our flow tube is a stainless-steel tube of diameter 6 cm and length 36 cm. 
The carrier gas was helium (at pressures near 1 Torr and temperatures near 295 K). 
Flow velocities ranged from 5000 cm/s to 6000 cm/s. Reagent gases are added to 
the flowing plasma through several movable and fixed inlets. A movable Langmuir 
probe 2 ® serves to measure electron densities rie as a function of axial position z. A 
quadrupole ion mass spectrometer analyzes the ion composition at the downstream 
end of the flow tube. 


Plasma generation and ion conversion 

The techniques of recombination measurements in flow tubes are quite well 
known and will not be discussed here. The method consists of creating a well- 
characterized plasma that contains only the desired ion species and measuring the 
rate at which free electrons disappear. The creation of the H3 + plasma involves 
several steps: 

(1) The carrier gas (helium) is partially excited to the metastable state and (to a 
lesser degree) ionized by passing it through a microwave discharge. 

(2) Argon (~ 15 % of helium) is added downstream to convert metastable helium 
atoms to Ar + ions by Penning ionization. After this, the plasma contains mainly 
Ar+ ions with a small admixture of He + and He2 + ions (< 0.1 % of Ar*). The 
electron density at this point is typically 4xl0 10 cm* 3 . 

(3) To create H3 + ions, H2 is added at densities from 10 14 to 10 15 cm* 3 . The 
formation of H3 + via the reactions 

Ar+ + H 2 Ar H + + H (+ 1.53eV) (2) 

ArH + +H2 — » Ar+ H3 + (+0.57eV) (3) 

(see Radzig and Smirnov 27 for thermochemical data) takes place very rapidly 
compared to recombination since the rate coefficients for reactions (2) and (3) are 
large (~ 10*9 cm 3 /s 28 ). 

^ Vibrational excitation and quenchfng of the H3 + ions 

If the entire energy of 2.1 eV that is released in reactions (2) and (3) were 
to appear as vibrational excitation, which seems unlikely, then H3 + could be pro- 
duced in vibrational states up to v=5. However, H3 4 - ions in v>l would be de- 
stroyed rapidly by proton transfer to Ar, 

H3 + + Ar — > Ar H + + H2 


( 4 ) 



which is energetically possible for H3 + ions with internal energies above 0.57 eV. 
The non-reacting, surviving H3 + vibrational states are the ground state, the V2=l 
bending-mode vibration at 0.3126 eV, and the vj=l breathing mode vibration at 
0.394 eV(see Lie and Frye 29 or Oka and Jagod 6 ). At [Ar]= 5x 10 35 cm' 3 , v>l 
ions would be destroyed by proton transfer (k ~10' 9 cm 3 /s) in less than 1 |isec, a 
time that is far shorter than the time scale of recombination (about 1 ms). 

The radiative lifetimes of the vj=l and V2=l levels (1.2 s and 4 ms, respec- 
tively) 19 are too long for radiative relaxation to be important Quenching of v=l 
ions in collisions with H2 is more effective. Amano's 333 absorption studies of the 
V2 hot band indicate that the V2=l level is quenched by H2 with a rate coefficient 
of approximately 3xlO' 333 cm 3 /s. This eliminates V2=l ions; they will be destroyed 
in 30 |is at [H2]= lxlO 14 cm' 3 . On the other hand, there is evidence 31 - 32 for 
slower quenching (k~ 1 O' 3 2 cm 3 /s) of other vibrational states. It seems very likely 
that the slowly-quenched ions are those in the vj=l level. A similar conclusion was 
reached by Bawendi et al. 33 who studied the absorption of H3 + hot bands in He/H2 
discharges. 

The foregoing considerations indicate that only H3 + (vi=l) ions live long 
enough to undergo recombination. Smith and Spanel 24 did not consider the de- 
struction of the higher vibrational levels by proton transfer with argon and conse- 
quently believed that most H3 + ions would be in states v>3. We do not agree with 
this conclusion; others have expressed the same criticism earlier 33 * 14 . This is an 
important point, since Smith and Spanel's proposed reconciliation of theory and 
experiment rests on the assumption that vibrationally excited ions dominate. 

There is no a priori reason to expect that para and ortho H3 + ions should 
recombine equally fast. According to Oka 34 , ortho-para equilibration by proton 
exchange with H2 is rapid. We assume that the proton-transfer rate is sufficiently 
large (lolO* 333 cm 3 /s would suffice at [H2] = 10 34 cm* 3 ) to keep the para-to-ortho 
ratio at a constant value in the recombining plasma. 

DATA ANALYSIS 


The most common method of deducing P from the slope of a graph of the 
reciprocal of the electron density ne(z) vs. axial position z (V being the plasma 
velocity). 


1 

n e (z) 


1 

n t (z = 0) 



( 6 ) 


provides quick, but only only approximate values of a z-averaged de-ionization 
coefficient. Since we are looking for deviations from a simple decay law, we were 
forced to employ a more rigorous approach (for details see ref. 25): It consisted of 
solving numerically the continuity equation for the electron density n e (r, z) in a 
flowing plasma with a parabolic velocity field. The computed solutions were then 
used to validate a computationally faster, approximate method that could be used 
for routine data analysis. The approximate method extracts a z-dependent de-ion- 
ization coefficient P(z) from the measured electron density n e (z). This form of 
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analysis was shown to be sufficiently accurate by using the exact model to gener- 
ate simulated n e (z) values and then analyzing the simulated data by the approxi- 
mate method. In addition, tests were made on experimental data obtained in 02 + 
plasmas. As expected, no systematic variation of (3(z) with z was found and the 
measured P agreed with the known recombination coefficient a(C>2 + ) . 

EXPERIMENTAL RESULTS 

40 decay curves n e (z) were recorded for [Ar]= 4x 10^ cm* 3 to 8x lO 1 ^ 
cm* 3 , and [H 2 ]=lx 10 14 cm' 3 to 1.7x 10 15 cm* 3 . Even the simplest form of data 
analysis (based on Eq. (9) ) showed that P increased noticeably with increasing H 2 
concentration (see Fig. 1). A similar, but less scattered graph was obtained when 
the average values of P(z) in the first 4 cm were plotted as a function of [H 2 ]. 

Figures 2 and 3 show examples of the electron density ne(z) and P(z) for 
two different H 2 concentrations. It is apparent that P is not constant but that it 
decreases with increasing z. This finding reproduces the observations of Smith and 
Spanel 24 ; the decay of the plasma clearly does not follow a simple recombination 
law. 

A smaller set of data was taken for D 3 + ions. The observations largely mir- 
rored those for H 3 + . Again, the recombination coefficient increased with increasing 
D 2 concentration, but the D 3 + de-ionization coefficients were significantly smaller 
(from 0.7xl0' 7 cm 3 /s to 1.3 xlO' 7 cm 3 /s in the range [D 2 ] = 1 to 14 x 10 14 cm' 3 ). 

INTERPRETATION OF RESULTS 

The interpretation should explain the enhancement of p upon hydrogen 
addition, the decline of P(z) with increasing position z, and if possible, it should 
elucidate the recombination mechanism. We attempted to explain the peculiar 
decay of the plasma by two different models. The first model assumes the presence 
of H 3 + ions in two different states. The second model invokes an entirely different 
explanation in terms of Rydberg molecules. 

Two - state model 

The assumption that the z dependence of P is caused by the presence of 
H 3 + ions in v=0 and in vi=l leads to 

P(z) = Oofo(z) + aifi(z) (6) 

where ao and cci denote the recombination coefficients of v=0 and vj=l ions, and 
foand fi are their fractional concentrations. Quenching of vj=l ions to v=0 by H 2 
(rate coefficient kq) and differing recombination losses make the concentrations z- 
dependent. It is not difficult to construct a kinetic model that yields a fit to the 
measured P(z) curve at [H 2 ] = 1.3 x 10 14 cm' 3 . Figure (2) shows fits for oto= 2.2 
xlO' 7 cm 3/ s, ai= 4 xlO' 8 cm 3 /s, f 0 = 0.7, and two quenching coefficients of 
kq=lx!0* 12 cm 3 /s and kq=lxl0' 13 cm 3 /s. However, when the same values are 



used to compute a fit to the data at [H2]= 1.7 x 10 15 cm* 3 ’ the model (see Fig. 3) 
indicates that (3(z) should now increase with z. A decline of (3(z) can be forced by 
reducing the quenching rate to kq=lxlO' 13 cm 3 /s, but the model fails to reproduce 
the observation that (3 at z=0 (only 1 cm from the reagent inlet)is larger at the 
larger H2 concentration. We eventually concluded that the two-state model does 
not provide a consistent explanation of the experimental findings. 

De-ionization by Rydberg molecule stabilization 

We now return to the subject that was mentioned in the Introduction, 
namely de-ionization mechanisms which involve interactions with third bodies. 
The |3(z) data suggest that the de-ionization coefficient is larger at higher electron 
densities and higher H2 densities. Is it possible that these particles affect the mech- 
anism of de-ionization directly, rather than changing the state of the recombining 
ions? At first sight, this does not seem to be viable proposition. Electron-stabilized 
and neutral-stabilized recombination are known to be too slow (by several orders 
of magnitude) 35 to have an effect on recombination in low-density (ne< 10 11 era' 3 ) 
plasmas. Third-body effects are inefficient because electron capture requires 
energy transfer to a third particle during the brief electron-ion encounter. 

In the following, we will argue that additional stabilization mechanisms 
exist in the case of molecular ions, especially those that do not undergo fast dis- 
sociative recombination. Theory indicates the H3 + ions may fall into this category. 
Non-dissociating, long-lived Rydberg states may be formed in collisions of 
electrons with molecular ions. These can be stabilized not only by energy- 
removing collisions but also by angular-momentum changing (/-mixing and 
possibly m-mixing) collisions. 

The following sequence may occur: In the first step, an electron of low 
angular momentum (s or possibly p) is captured into a ro-vibrational level of a 
Rydberg state of H3, 

e- + H 3 + H 3 ** (v=l,n,s). (7) 

The n=7 state should be favored since excitation to the V2=l state (0.3126 eV) in 
n= 7, 1=0 is resonant at an electron energy of +0.0309 eV (1.2 kT at 300 K). 

We postulate that in the second step the angular momentum of the Rydberg 
molecule is randomized in collisions with ambient electrons (1-mixing and possibly 
m-mixing). It is known from studies of dielectronic recombination 36 and zero- 
electron-kinetic-energy (ZEKE) spectroscopy 37 that redistribution of angular mo- 
menta of Rydberg electrons (1-mixing) due to electric microfields and collisions 
can increase predissociation and autoionization lifetimes by several orders of mag- 
nitude. The basic reason is that the wave functions of high-1 states do not interact 
with the ion core. Since the lifetimes of some H3** autoionizing states are quite 
long (of the order of 10' 7 s) 38 > 39 and the 1-mixing efficiency of thermal electrons 
is extremely large (rate coefficients >10' 3 cm 3 /s) 4 ®, 1-mixing should occur in less 
than 10' 7 s under the conditions of our experiment (n e = 4 xlO 10 cm' 3 ). 
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In the third step, the long-lived H 3 **(n=7, v=l, 1>I ) can be returned to the 
autoionizing state (s, or p) by further collisions, they may be quenched to lower n- 
states by neutral atoms or electrons, they may undergo slow predissociation, 
radiate to a lower state, or react with an H2 molecule. Some of the decay rates are 
known or can be estimated. Collisional energy loss (n-changing) with helium 
atoms is not important 41 compared to radiative decay. The rate coefficient for n- 
changing collisions of electrons with n=7 Rydberg atoms (electron quenching) has 
been measured by Devos et al. 43 to be -10*5 cm 3 /s. Radiative decay rates for n=7 
hydrogenic states are close to lxlO 6 s* 1 for s-states , 2x 10 7 s* 1 for p-states, and 
about 2x10^ s* 1 for a statistical mixture of 1-values up to 1=6 42 . The reaction with 
H 2 

H 3 **+H 2 ->H 2 + H 3 * (8) 

provides the needed enhancement of the loss rate of electrons in the presence of 
H 2 . It will compete with radiative decay if it occurs on a time scale of <l|is . At 
[H 2 ] = lx 10 15 cm* 3 , this requires a rate coefficient of - lxl0* 9 cm 3 /s for reaction 
(8), which is not unreasonable. 

The set of reaction equations is easily solved to obtain an effective de- 


ionization coefficient. The capture and autoionization step 

H 3 + + e * <-» H 3 ** (n=7, s,p) (9) 

is described by a capture coefficient kc and an autoionization rate v a . Stabilization 
by 1-mixing due to electrons 

H 3 ** (n=7, s, p)+ e * H 3 ** (n=7, l>l)+e * (10) 

and its reverse are described by rate coefficients k s and ki, respectively. The fol- 
lowing destruction reactions for H 3 ** were included: 

H 3 **(n=7) -» H 3 *(n<7) + hv (radiation, rate Vrad) (Ha) 

H 3 ** H 2 +H (or 3H) (predissociation, rate v p ) (lib) 

H 3 **(n=7) + e* — » H 3 *(n<7) +e* (rate coefficient kqe) (11c) 

113 **+ H 2 — » H 2 4 -H 3 (lown) (rate coefficient kn ) • (lid) 


To allow for the possibility that the rate coefficients for H 3 ** in the 
autoionizing states (s, p) differ from those in higher 1-states, we treat their de- 
struction separately and attach the subscripts "0" and "1" to the rate coefficients. 
After some algebra, the effective de-ionization coefficient can be put in the form 


where 

and 


a _ 1 dn t _ 

*"l + 

(nA +V t )\r a V a 

Vo = v ra( j,o + Vp t o + kqe ne + kn [ H 2 ] 
v l = Vrad i + Vp t i+ kqe n e + kn [ H 2 ]. 


( 12 ) 



It is possible to find a set of plausible rate coefficients so that Eq. (12) 
gives a reasonable fit to the data. Fig. 4 shows a comparison to experimental val- 
ues of p (n e ) for the following set of rate coefficients: 

k c = 4.4x 10 ' 7 cm 3 /s, v a = 1.3xl0 8 s _1 , k s = 1.25 xlO " 2 cm 3 /s 
kj =3.1 xl(H cm 3 /s, Vra^o =Vrad.l = 1-1 xl^S ' 1 
Vp,o=Vp,i =0, kn = 4.5 x 10 ' 9 cm 3 /s, kq e = lx 10 ’ 5 cm 3 /s 
We do not claim that this set of model parameters is unique or that the 
agreement with the data is sufficient proof for the validity of the model. More de- 
tailed studies would be required to prove that this model is valid. It is interesting, 
however, that the fitting parameters are quite close to those that were estimated 
earlier. If this model is correct, it would solve the theoretical dilemma of having 
to account for fast "binary recombination" of electrons with H 3 + ions since this 
process would not be needed to explain the observed plasma de-ionization. 

Our model is compatible with the most recent merged-beam observations 
of Mitchell and coworkers 43 . They found that the apparent H 3 + recombination 
cross sections decreased by about a factor of five when a "deflection field” in the 
post-collision region was increased in magnitude. The effect is ascribed to field-in- 
duced re-ionization of autoionizing H 3 ** molecules that were formed in electron- 
H 3 + collisions. The role of electric fields in such experiments is not yet perfectly 
understood, but the observations demonstrate that long-lived H 3 ** states can be 
produced in such collisions. The nature of these states and the reasons for their 
unexpectedly long life-time, however, arc not at all clear. 

CONCLUSIONS 

We reach the tentative conclusion that the recombination of H 3 + plasmas 
does not exclusively occur by binary recombination. If this is true, then the results 
of afterglow measurements may be compatible with the theoretical expectation that 
binary recombination of H 3 + ions in v =0 should be rather slow. 

The "true" binary recombination coefficient of H 3 + is needed for some 
applications, e.g. for models of the interstellar medium. It is conceivable that the 
recombination cross sections measured in merged-beam or ion-storage ring exper- 
iments are close to the true values, but the effects of residual electric Fields in such 
experiments need to be examined more closely. 
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